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Abstract 

The time evolution of vector-meson spectral-functions is studied within a kinetic 
theory approach. We implement this formalism in a BUU type transport model. 
Applications focus on p and uj mesons being important pieces for the interpretation 
of the di-electron invariant mass spectrum measured by the HADES collaboration 
for the reaction C + C at 2 AGeV bombarding energy. Since the evolution of the 
spectral functions is driven by the local density, the in-medium modifications are 
tiny for small collision systems within this approach. 

1 Introduction 

Di-electrons serve as direct probes of dense and hot nuclear matter stages during the 
course of heavy-ion collisions [Tj [2] . The superposition of various sources, however, re- 
quires a deconvolution of the spectra by means of models. Of essential interest are the 
contributions of the light vector mesons p and uj. The spectral functions of both mesons 
are expected to be modified in a strongly interacting environment in accordance with 
chiral dynamics, QCD sum rules etc. [H EJ HI [5] . After the first pioneering experiments 
with the Dilepton Spectrometer DLS [6] now improved measurements with the High- 
Acceptance Di-Electron Spectrometer HADES [3, EJ [9] start to explore systematically 
the baryon-dense region accessible in fixed-target heavy-ion experiments at beam ener- 
gies in the 1-2 AGeV region. The invariant mass spectra of di-electrons for the reaction 
C + C at 1 and 2 AGeV are now available [6j [8j [9] allowing us to hunt for interesting 
many-body effects. 



There are several approaches for describing the emission of real and virtual photons 
off excited nuclear matter: 

(i) A piece of matter in thermal equilibrium at temperature T emits e + e~ pairs with 
total momentum q at a rate dN/d A xd 4 q = a 2 / (M 2 7r 3 )fBlinIl em , where fs is the bosonic 
thermal distribution function and H em denotes the electromagnetic current- current cor- 
relator, Il em = —% J d 4 xe iqx ((Tj f *(x)j fl (0))), the imaginary part of which determines 
directly the thermal emission rate. Here, a stands for the electromagnetic fine structure 
constant, M means the invariant mass of the di-electron, T is the operator for time 
ordering which acts on the electromagnetic current operator j^. This rate may be com- 
bined with a global dynamic model which provides us the space averaged quantities as a 
function of time [10] or even adopting also a time average [TT] . 

(ii) Some sophistication can be achieved by employing a detailed model for the space- 
time evolution of baryon density and temperature, e.g., as delivered by hydrodynamics. 
One may also extract from transport models such parameters, where, however, local off- 
equilibrium and/or anisotropic momentum distributions hamper a reliable definition of 
density and temperature. Nevertheless, once the rate is given one has a very concise 
approach, as realized, e.g., in [12]. A similar approach has been presented in [T3] . 

(iii) Microscopic or kinetic transport models do not require isotropic momentum distri- 
butions or local equilibrium. Once general principles are implemented, transport models 
also provide a detailed treatment of the emission of electromagnetic radiation in heavy-ion 
collisions. 

Our approach belongs to item (iii). The time evolution of single particle distribution 
functions of various hadrons are evaluated within the framework of a kinetic theory. 
We focus on the vector mesons p and u. The p meson is already a broad resonance 
in vacuum, while the u meson may acquire a noticeable width in nuclear matter |14j . 
Therefore, we are forced to treat dynamically these resonances and their decays into 
di-electrons. Resorting to consider only pole-mass dynamics is clearly insufficient in a 
microscopic approach [T5]. Instead, one has to propagate properly the spectral functions 
of the p and u mesons. This is the main goal of our paper. We consider our work as being 
on an explorative level, not yet as a firm and deep theoretically founded prescription of 
dealing with di-electron emission from excitations with quantum numbers of p and u 
mesons off excited nuclear matter. 

Our paper is organized as follows. Essential features of our transport model are 
outlined in section 2. In Subsection 2.1 we describe how the mean field potentials enter 
in the relativistic transport equation. Subsection 2.2 introduces the dynamics of broad 
resonances and its implementation in the test-particle method. The crucial quantities 
for the spectral functions are the self-energies dealt with in subsection 2.3. Subsections 
2.4 and 2.5 are devoted to particle production and di-electron emission, respectively. 
Numerical results of our simulations and a tentative comparison with published HADES 
data are presented in sections 3 (2 AGeV) and 4 (1 AGeV). Discussion and summary can 
be found in section 5. 

The present analysis supersedes [16] - Further analyses of HADES data have been 
performed in [13 O [13 ED]. 
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2 Treatment of heavy-ion collisions 



2.1 The standard BUU treatment 

The employed BRoBUU computer code for heavy-ion collisions developed by a Budapest- 
Rossendorf cooperation solves a set of coupled Boltzmann-Uhling-Uhlenbeck (BUU) equa- 
tions in the quasi-particle limit |21j 



dFi dHdFi dHdF t ~ r- — — -— - 

aT+aF&F- ^^ = I>> ^=v / K + ^(p,x)) 2 + p 2 (i) 

for the one-body distribution functions i^(x, p,t) of the various hadron species i, each 
with rest mass m», in a momentum and density dependent mean field U. The scalar mean 
field U is chosen in such a manner that the Hamiltonian H equals H = \J mf + p 2 + U™ r 
with a potential [7™ r calculated in the local rest frame as 

IP" - A n j.r(" Wr 2 f ^ f "(x,p') m 
U ' - A n +B {nJ +C nJ J2^ 1 + (2 ^f (2) 

where the parameters A, B, C, r, A define special types of potentials, while n, n and 
Fn stand for the baryon number density, saturation density and nucleon distribution 
function. We use the momentum dependent soft potential defined by v4=-0.120 GeV, 
5=0.151 GeV, r=1.23, C=-0.0645 GeV, A=2.167 GeV. The BRoBUU code propagates 
in the baryon sector the nucleons and 24 A and N* resonances and additionally tt, n, a, uj 
and p mesons. Different particle species are coupled by the collision integral Cy which also 
contains the Uhling-Uhlenbeck terms responsible for Pauli blocking of spin-1/2 hadrons 
in the collision as well as particle creation and annihilation processes. 

The set of coupled BUU equations is solved by using the parallel-ensemble test-particle 
method [22, 23J, where we introduce a number of parallel ensembles. In each ensemble a 
test particle represents a real particle (nucleon, resonance, pion etc.); collisions happen 
only within the same ensemble. On the other hand, when calculating such quantities 
as densities, Pauli blocking factors etc. we average over the ensembles in each time step. 
This method transforms the partial differential-integro equations (pQ) into a set of ordinary 
differential equations (looking like equations of motion) for a number of test particles. A 
default version of the code has been applied to strangeness dynamics [24] . 



2.2 Off-shell transport of broad resonances 

Recently theoretical progress has been made in describing the in-medium properties of 
particles starting from the Kadanoff-Baym equations [25] for the Green functions of 
particles. Applying first-order gradient expansion after a Wigner transformation one 
arrives at a transport equation for the retarded Green function [26], |27] . In the medium, 
particles acquire a self-energy H(x,p) which depends on position and momentum as well 
as the local properties of the surrounding medium. They have a finite life time which is 
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described by the width T related to the imaginary part of the self-energy. Their properties 
are described by the spectral function being the imaginary part of the retarded propagator 
A(p) = —2ImG ret (x,p). For bosons the spectral function is related to the self-energy via 

A(p) = ^ P , (3) 

{E 2 -p 2 - m 2 - ReZ ret (x,p)) 2 + \T{x,p) 2 

where the resonance widths V and f obey f (x,p) = —2ImY7 et « 2moT, and mo is the 
vacuum pole mass of the respective particle. The spectral function ([3]) is normalized as 
J dp 2 A = 2n. (We omit here the label denoting the particle type for simplicity reasons.) 

To solve numerically the Kadanoff-Baym equations one may exploit the above test- 
particle ansatz for a modified retarded Green function (see Refs. [26| 127]). This function 
can be interpreted as a product of particle number density multiplied with the spectral 
function A. The spectral function can significantly change in the course of the heavy-ion 
collision process. Therefore, the standard test-particle method, where the test-particle 
mass is a constant of motion, must be extended by treating the energy E = p° of the 
four-momentum p as an independent variable. 

Equations of motion for test particles follow from the transport equation. We use the 
relativistic version of the equations which have been derived in Ref. [26J: 

"* 1 1 f 2f + O pR e^ + ' n2 -<- ReE "'§ r t], (4) 



dt i - c V. r 

dp - 11 (d x Re^ + m2 -<7 ReJ:ret d x r), (5) 



dt 1 - C 2£ V r 

dE 1 1 W-'+" ,2 -" 1 °r Jtesr Vi, ( 6) 



dt 1 - C 2E v r 
with the renormalization factor reads 

C= ±(s ERe ^> + m »- m \- ReEr * 8 E fy (7) 

In the above, m = y ' E 2 — p 2 is the mass of an individual test-particle number of a given 
hadron specie. The self-energy E re * is considered to be a function of density n, energy E, 
and momentum p; thus the dependence on time and position comes only from its density 
dependence. Partial derivatives with respect to any of the four variables, t,x,E,p, are 
understood taking the three other ones fixed. The quantity 1 — C has been introduced 
to ensure that the test particles describe a conserved quantity [27] . 

The change of the test-particle mass m can be more clearly seen combining Eqs. ([51) 
and © to 

dm 2 1 ( d R ^ ret | m 2 - m 2 - ReL ret d A ^ 



dt 1 - C \dt r dt 

with the comoving derivative d/dt = dt +p/Ed x . This equation means that the square 
of the particle mass tends to reach a value shifted by the real part of the self-energy 
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within a range of the value of f . Thus, the vacuum spectral function is recovered when 
the particle leaves the medium. This ensures the smooth transition from the in-medium 
behavior to the vacuum properties. 

The equation of motions of the test particles have to be supplemented by a collision 
term which couples the equations for the different particle species. It can be shown [27] 
that this collision term has the same form as in the standard BUU treatment. 

The off-shell transport has been implemented in simulations for the propagation of 
p and us mesons to study the di-electron productions in and pA reactions. Early 
approaches E2] did not automatically provide the correct asymptotic behavior of the 
spectral function and an auxiliary potential were introduced to cure this problem. The 
above equations of motion, which do not have this deficit, were applied to di-electron 
production in pA collisions in [30] and for AA collisions in [T7J . 



2.3 Self-energies 

To solve the Eqs. (jUED one needs the knowledge of the self-energies. Here one faces the 
need to decide which effects to take into account in the expression for the retarded self- 
energy Y7 et in the medium. That is because the BUU transport equations themselves 
already contain some part of in-medium effects that usually are considered in theoretical 
models in local density and local equilibrium approximation [2H E21 E3]- For instance 
models for in-medium effects of p mesons usually take into account the iV(1520)-nucleon- 
hole loop for the self-energy, the corresponding vertices are accounted for in BUU via 
p-nucleon scattering and absorption through the iV(1520) resonance. 

In our calculations we employ a simple schematic form of the self-energy of a vector 
meson V: 

Ti 

ReSL e * = 2m v Am v — , (9) 

n 

lmE? = my(rr + ^==). (10) 

V 1 — v z 



Equation ([9]) causes a "mass shift" Am = y/m v + ReS^f — my characterized by Amy 
and being roughly proportionally to the density n of the surrounding matter. The imag- 
inary part contains the vacuum width T™ c the energy dependence of which is described 
by a form factor [3l] . The second term in Eq. (fTUj) results from the collision broadening 
which depends on density, relative velocity v and the cross section oy of the vector meson 
in matter. This cross section ay is calculated via the Breit-Wigner formula 

= 4tt v 2 J R + 1 5iv, fl rff 

av dv 3 ( 2J ' +1 )( s - m «) 2+s ( r ") 2 

for forming resonances with masses tor, angular momenta Jr, partial widths Ty^, total 
widths r^* with energy y/s and relative momentum q in in the entrance channel. In 
vacuum the baryon density n vanishes and the resulting spectral function A vac is solely 
determined by the energy dependent width Ty 20 . We remark that the decay of a test 
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particle is determined by its vacuum width. The life time of a test particle is, furthermore, 
reduced by the absorption during the two-particle collisions which characterizes the total 
width ImS^. 

In our actual numerical implementation we assume that the spectral function of the p 
(u) meson vanishes below two (three) times the pion mass, respectively (For a discussion 
of this issue see, for instance, [35] •) If a p meson is generated at normal nuclear matter 
density n its mass is distributed in accordance with the spectral function (see Eq. ffT6l) 
below). If the meson propagates into a region of higher density then the mass will be 
lowered according to the action of ReS re * in Eq. (jHJ). However if the meson comes near 
the threshold the width T becomes small and the second term of the right hand side of 
Eq. (jHJ) dominates and reverses this trend leading to an increase of the mass. 

The life time of unstable particles is also accessible in the framework of the transport 
equations for resonances. As it was shown in Ref. [36J this description leads to a life time 
r = d5/dE, where 8 is the energy dependent scattering phase in the formation or the 
respective decay of the resonance. Although this relation is known for a long time [37] 
it was introduced only recently in the context of a BUU transport treatment |38j . This 
prescription is very different from the commonly used formula r = h/T. Especially if the 
resonance is a p wave resonance the life time tends to small values near the threshold 
in the former case, while it approaches large values in the latter one. If the resonance 
decays into several channels, the total width is the relevant quantity which describes the 
phase of the amplitude common to all decay channels: 

-if 

tan 8 = — tt 1 =— r 

p 2 - ml- ReS re * 

Therefore, the decay rate into a special channel c is given by the partial width T c = b c r tot 
according to 

T -i = bcT tot = bc ( dS / dE yi ( 13 ) 

with b c being the branching ratio of the decay into channel c. If we do not mention 
otherwise, we use the standard prescription for the life time, but in some cases we study 
the effect of using Eq. (jT3"]l . as well. 

2.4 Particle production 

In most instances, a vector meson V is created by the decay of a baryon resonance R 
in the BRoBUU code. Thus, mesons are created in two-step processes like NN «-> NR 
with subsequent decay R «-> VN. As mentioned above the BRoBUU model includes 24 
non-strange baryon resonances. Their parameters (mass, width and branching ratios) are 
determined by a global fit to pion-nucleon scattering data, while resonance production 
cross sections are fitted to inelastic nucleon-nucleon scattering cross sections [31] • Since 
there are very few np —>■ RN data we assume that 

&np^RN — dCpp^RN, 

(14) 

6 



(12) 



where a is a channel independent (except one, see below) constant, and its value a = 1.34 
is obtained from a fit to the few existing data. We use that prescription for all resonances 
(R) except for N(1535), the main source of rj meson, where experimental data indicate a 
much higher value of 0^(1535) = 5. 

Our approach is in contrast to other ones where individual elementary hadron reaction 
channels are parameterized independently from one another. Using such coupled channel 
approach could allow us to obtain cross section for not or poorly measured channels. 

The in-medium spectral functions of u and p mesons also have to be taken into account 
when their test particles are created. In the resonance decay the mass distribution of the 
generated test particles for mesons results from an interplay of phase-space effects and 
the in-medium spectral functions A of the created meson. For the decay of a resonance 
of mass rriR in a meson of mass m and a baryon of mass we use the phase space 
distribution in the final state with a constant matrix element squared \A4\ 2 



T = Af d A p N 8(p 



m 



d 4 Pv —A(pv)\Mf 
In 



(15) 



from which the distribution 

dN R->NV 



dm 



v 



Af m v \ l l 2 {m 2 R) m 2 N , my) A(m v ) 



(16) 



results, where A is the triangle function A(a 2 , b 2 , c 2 ) = (c 2 — a 2 — b 2 ) 2 — (2ab) 2 . Af is an 
appropriate normalization factor. 

We also include meson emission during a transition R — > R'V from a resonance state 
R to another resonance R! with R' = A(1232), JV(1440), iV(1520), JV(1535). 
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Figure 1: Production cross sections of u (left), p° (middle) and rj (right) mesons in pp 
collisions as a function of the excess energy in comparison with data fSQ \41\ \42 
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In Fig. [T] we compare the cross section calculated with our parameter set with data 
measured by the collaborations SATURN, COSY and DISTO [33 EH HD H21 H31 Hfl 
1351 1361 137] . The relevant range region for collisions in the 1-2 AGeV region is at 
excess energies below 0.5 GeV. We recognize that the u production in pp collisions is well 
reproduced by our model parameters. The pn — ► pnui cross sections are about 1.5 times 
larger than the pp cross sections. The one-boson exchange model in [4_8] predict even a 
ratio of two. 

For the p° production near threshold there are not many pp measurements near thresh- 
old. Measurements are hampered by the large p width which make it difficult to discrim- 
inate p mesons from sequential two pion emission. At an excess energy of 0.33 GeV a 
p° cross section of 23±9 pb [13] has been measured, where p° mesons were identified by 
pion pairs with with masses above 0.6 GeV. The here employed global fit including many 
elementary channels overestimates this cross section by a factor of 2.3. As a consequence 
we will reduce the p production by this factor in the following. 

With respect of n production our model describes well the production in pn collision, 
however seems to underestimate the production in pp collisions. Since the cross sections 
in pp collisions are anyhow smaller than those of pn collisions, this fact will not seriously 
affect our final results. 

Furthermore the p mesons can also be created in pion annihilation processes tt + tt — > p 
(see below). 



2.5 Di-electron production 

The di-electron production from direct vector meson decays V — > e + e~ is calculated 
by integrating the local decay probabilities along their trajectories in accordance with 
Eq. (1131) . The branching ratios b c of the vector mesons are taken from experimental data 
at their pole masses. The mass dependence of this branching ratio is assumed to behave 
proportional to m^ 3 in accordance with the vector meson dominance model. 

The subleading so-called direct channel tttt — > p — > e + e~ is treated with the p meson 
formation cross section 

a(M) = ^2m° p V(p)A p (17) 

with m p being the pole mass of the p meson and T p (p) the vacuum width of the p 
resonance. The in-medium effects are encoded in the spectral function A. 
The p meson produces di-electrons with a rate of 

,o\ 3 

b c T{p). (18) 



dN P+P - ( m p 



dt 

We also include into our simulations a bremsstrahlung contribution which is guided 
by a one-boson exchange model adjusted to pp virtual bremsstrahlung and transferred to 
pn virtual bremsstrahlung [49] . Actually, we use 

da <j j_ a 2 f d 3 qR 2 (s) 



dM M 6vr 3 J g 3 R 2 (s 
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Here M is again the e + e~ invariant mass, R2 denotes the two-particle phase space volume, 
\fs stands for the c.m.s. energy, s is the reduced energy squared after the di-electron 
emission, and o~±(s) is the transverse cross section. Equation ffl9|) can be approximated: 



da 



or o tot s - (mi + m 2 ) 



dM 3tt 2 M 



ln( 



Qrnax QOri 



-) 



qor, 



(20) 



This approximation is applied to pn and 7rN collisions using the respective corresponding 
total cross sections o tot \ e cm stands for the energy of the charged particle in the rest system 
of the colliding particles with masses mi and m 2 , qomax = (s + M 2 — (mi + m 2 ) 2 )/2y / i 
is the maximum di-electron energy, and q max = \J q\ max — M 2 denotes the maximum 
di-electron momentum. It should be noted, however, that this cross section is still rather 
uncertain. 

An essential di-electron contribution comes from the Dalitz decays of ir°, rj, uo mesons 
and the excited baryon resonances emitting a di-electron together with a photon or nu- 
cleon. The decay rate Yoai for a di-electron of mass M for mesons can quite generally 
brought into the form 



dT 



Dal 



dM 



_Aa_ 

3ttM 7 



1 - 



M 5 



mf 



F(Mf 



(21) 



The Dalitz decay rates are assumed to be given by the photon partial T 7 width which 
have been taken from experiment [50J. The relevant form factors F(M) for the mesons 
being considered are summarized in [51]. 
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Figure 2: Dalitz- decay of A(1232) for two energies (masses) m&: 1.232 GeV and 1.85 
GeV with the description from above 123] \52\ [55j [37]/ (left panel) and the spectral function 
with two different cut-off' prescriptions (right panel). See text for further details. 

The Dalitz decay of the baryon resonances is treated as in [52]. The most important 
contribution to the di-electron spectra of these come from the A(1232) resonance. There 
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are also other models [511 [231 [53] for the Dalitz-decay. As can be seen on the left panel of 
Fig. [21 these models agree very well for resonance decays from the peak mass, however 
they differ substantially for A resonances with energies (masses) relevant for studying the 
vector meson region. There is another uncertainty concerning these high-energy (mass) 
A(1232) resonances. The width and consequently the spectrum of these resonances are 
sensitive on the cut-off for high masses. Here we show two possible parameterizations: 
one from Moniz [53] and the other one from Manley used in the Particle Data Book |50j . 

The number of A(1232) resonances at energy (mass) around 1.85 GeV may depend 
on the cut-off prescription by a factor of 3. Their Dalitz decay (see left of Fig. [2]) 
panel may differ by a factor of 4. So the Dalitz decay contribution of the A(1232) 
resonance is uncertain by more than an order of magnitude in the vector meson region. 
This uncertainty may only be clarified by a detailed comparison of the calculation for 
pp — > ppe + e~ with forthcoming experimental data. Using their angular dependence one 
can localize the different channels and then fix their magnitude. Here we would like to 
mention that different groups use different prescription for that channel which it is one 
reason why the predictions, especially for the A(1232) Dalitz decay contribution, are 
different. 



3 Results for 2 AGeV 



3.1 Spectral function dynamics 

We employ the above described code for the reaction C(2 AGeV) + C, where data from 
HADES are at our disposal [8J. In the present explorative study we are going to contrast 
simulations with and without medium modifications of p and u mesons to elucidate 
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Figure 3: Influence of medium effects on the spectral functions for p (left panel) and uj 
(right panel) mesons at various nuclear matter densities and velocities in matter. 
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to which degree medium effects may become visible in the light collision system under 
consideration. In doing so we use fairly schematic medium effects (to be considered as an 
upper limit) condensed in a "mass shift" described by the above parameter Am^ = —50 
MeV in Eq. (J2J) for the u meson. Previous CB-TAPS data [55] suggested indeed such a 
mass shift. (See, however, [56J for a critical discussion of this data.) This problem 
is also investigated experimentally [57] and theoretically [31] [58] [59] in calculating the uj 
spectral function. The use of QCD sum rules [60] then can be utilized to translate this 
shift into a significantly larger shift for the p meson (dictated essentially by the Landau 
damping term); we use here Am p = —100 MeV. We are aware of experiments as reported 
in Ref. [61] which do not observe a noticeable shift of the p meson excitation strength. 
Nevertheless, several theoretical attempts are made to predict a possible p "mass shift" 
during the last decade. Many of them predict a fairly large shift of strength of p excitation 
to lower energy [3] H] , see also [32] [33] [35] [62] [63] Ell [65] . Thus, we keep this (presumably 
too large a) value to illustrate whether it would have a significant imprint on the observed 
spectra. 

The spectral function for p and u mesons are shown in Fig. [3] at two different densities 
of nuclear matter and two meson velocities in comparison with the vacuum spectral 
function. We would like to emphasize the strong velocity dependence of the widths, in 
particular for the u mesons. Despite the excitations of various baryon resonances the 
spectral functions appear as relatively smooth distributions. 
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Figure 4: Time evolution of the masses (upper panels) of about 100 test particles of uo 
mesons in a C + C collision at 2 AG 'eV 'kinetic beam energy at an impact parameter of 1 
fm. The red open circles indicate the time instant when test particles are annihilated due 
to either pionic decay or resonance absorption. The lower part displays the corresponding 
local densities for the same test particles. The left (without mass shift) and right (with 
mass shift) panels indicate the effect of the spectral function and decay properties. 
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Let us now consider the effect of the mass evolution of the u and p vector mesons 
given by the equations of motion Eqs. (jUl6]). The ensemble of test particles of mesons is 
generated in dense matter where their masses are distributed in accordance with their 
broadened and mass shifted spectral function. In Fig. H] we show the time evolution of 
a small ensemble of u test particles. In the calculations the nuclei touch each other at 
a time of 2.5 fm/c while the density peaks at about 6 fm/c and drops at 8 fm/c below 
saturation density (see lower part of the figure). At maximum density most of the vector 
mesons are created, afterwards the mass distribution gets narrower. Only a few of the 
u mesons decay in the dense phase where their masses deviate strongly from the pole 
mass value. If a low density is reached the vacuum spectral function dominates the di- 
electron decays which leads to a sharp peak at the pole mass. The right hand part of 
Fig. H] shows the test particles with a spectral function with both the mentioned mass 
shift and collision broadening, while the distribution on the left hand part is calculated 
with collision broadening only. 
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Figure 5: Time evolution of the masses of about 400 test particles of p mesons in a 
central C + C collision at 2 AGeV kinetic beam energy at an impact parameter of 1 fm. 
The particles in the left panel are calculated with a spectral function without a mass shift 
while the middle and the right ones the indicated mass shift were used. Note the different 
prescription for the life time indicated in the legends. 



Figure [5] displays the analog behavior of the test particles of p mesons. In the high 
density stage one recognizes the large spread of the p masses of about 300 MeV due to 
the imaginary part of the self-energy. Most of the p mesons decay rapidly (life time r < 2 
fm/c). If we use the relation r = 1/T for the life time then essentially the high-mass 
particles with their large width decay rapidly and, hence, have little chances to radiate 
di-electrons during their short life time. They look like flashes occurring only in a narrow 
time interval, thus not causing a longer paths of adjacent points in the mass vs. time 
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plot. New mesons are readily created at later times and lower densities. A few low- mass 
p mesons survive these periods and can still be found at 25 fm/c. Therefore, one expects 
a shift of the di-electron spectra to lower masses. Quite a different picture shows the 
right hand panel where the life time is calculated accordingly to Eq. ffT31) . Here the low 
mass particles have a shorter life time than the more massive ones. 

In Fig. [H] we exhibit the resulting di-electron spectra from p and uj decays for the 
different prescriptions of the life time. The dashed line shows the p and uj spectra if 
the masses of the test particles are kept fixed to the values at the instant of creation 
(green dash-dotted line). In this case the di-electron spectrum reflects the initial mass 
distribution which contains the high density spectral function. The effect ist best visible 
for uj mesons near the pole mass: here, the peak at the pole mass would nearly disappear 
when disregarding the mass evolution. The time evolution of the off-shell propagation 
pushes the resulting di-electron spectra towards their vacuum spectral function. If the 
life time of the vector mesons follows the standard expression r = K/Y then the low-mass 
vector mesons have sufficient time to reach their pole mass. This behavior is also clearly 
seen (solid line) for the uj meson. However, if Eq. ffTB"]) controls the decay this shift is 
hindered by the earlier decay of the low mass mesons (compare solid and dotted lines). 
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Figure 6: Di-electron spectra from direct decay of uj (left panel) and p (right panel) 
mesons calculated with different assumptions for the dynamics and the spectral functions, 
"vacuum" : vacuum spectral function, "matter" : collision broadening and mass shifts with 
two differen assumption of the life time, "static": in-medium spectral function while the 
mass evolution of the mesons is switched off. 



Finally we investigate the distribution of the emitted di-electrons as function of the 
density of the emitting region. We consider the effect of the uj mesons in three different 
density regions: (i) the density n < n /3, (ii) n /3 < n < n , and (iii) n > n . For the 
light CC system di-electrons from all region have similar masses and can therefore hardly 
disentangled in experiment, see Fig. [3 
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Figure 7: Contribution to the di- electron yield from to mesons in various density regions 
compared with the total yield (thick solid line labelled by "total"). The left picture is 
calculated with the standard life time for u mesons, while the right panel shows the effect 
when using the life time r = dbjdE. 



3.2 Comparison with HADES data 

While Figs. HI [5] and [6] refer to the central point of our work, we now look how the evolving 
spectral functions compare with data. In doing so the other di-electron sources have to 
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Figure 8: Di-electron invariant mass spectra for C(2 AGeV) + C calculated with in- 
medium spectral function. Individual contributions are depicted. Left panel: Full phase 
space. Right panel: With experimental filter IbUtf and compared to HADES data J2J/. 
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Figure 9: As Fig. but with vacuum spectral function (collision broadening is included). 
The contribution from uj and p mesons are indicated. The violett dashed line shows the 
spectra calculated with the in-medium spectral function from Fig. [3 

be included. The obtained di-electron spectra are represented in Figs. [8] and [9j In Fig. [8] 
the results are exhibited obtained by including the above described pole-mass shifts as an 
additional medium modification of p and uj mesons. In Fig.[9]only the collision broadening 
is employed which modifies the imaginary part of the spectral function. For comparison 
with the data, the HADES filter has been applied [66] accounting for the geometrical 
acceptance, momentum cuts and pair kinematics. The filter causes a reduction of the 
strength and a smearing of the invariant masses of the di-electrons. The result of this 
filtering is always shown on the right hand panel of the figures. 

In these figures we show various contributions to the di-electron rate. Important 
low-mass di-electron sources are ir° and 77 Dalitz decays which are proportional to the 
multiplicities of their parents. The TAPS collaboration has measured [67] the ir° and rj 
production cross sections of 707±72 mb and 25±4 mb which have to be compared to our 
calculations of 870 mb and 23 mb in the same reaction at the same energy. While the 
values for pion production are overestimated the 77 production is quite in agreement with 
the data. (Note that the presently employed cross sections rely on a global fit of many 
elementary reactions which is not optimized for special channel.) The Dalitz decays of p 
and uj mesons and nucleon resonances do not contribute noticeably. 

Comparing Figs. [8]and[9j the mass shifts of the vector mesons do not have a noticeable 
effect on the overall shape of the di-electron spectra although the peak position of the 
p mesons is clearly shifted from 0.7 to 0.5 GeV. However, the large contribution of the 
cocktail of the other sources cover the effect of the p mesons. Furthermore most of the 
uj mesons decay outside the dense zone and are therefore not very sensitive for medium 
effects. Since the fine structure {uj peak) is not yet resolved in the data a conclusive 
decision cannot be made. 
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Progress could be made if the di-electron mass resolution is improved to identify the 
uj peak. However, our calculations do not point to the possibility of a two-peak structure 
(resulting from a superposition of vacuum decays and in-medium decays) or a substantial 
smearing of the to peak due to a density dependent shift analog to the consideration of 
the cf) meson (see [SE])- 

The present set-up provides a reasonable description of the HADES data [7] for di- 
electron masses below 0.6 GeV. In the higher mass region some overestimation of the 
data is recognized. With respect to the uncertainties of the cross section pp — > ppp at 
threshold and generally the pn channel as well as the role of the resonance channels one 
we could try to improve the agreement by rescaling the p and uj contributions. In doing 
so we assume that the spectral shapes are unaltered. To get a better agreement with 
data one needs to decrease artificially the cross sections for vector meson production by 
factors 0.2 for uj and 0.8 for p. Figure [TU1 exhibits this "optimized" comparison with data. 

The transverse momentum spectra for three invariant mass bins are exhibited in 
Fig. [TTJ One recognizes a good agreement with these mult i- differential data. 



3.3 Effects in larger collision systems 



The C + C system is rather light and a consequence of it is that the maximum density 
is about 2.5 n (see Fig. UJ. A heavy system has a longer living high-density stage 
reaching densities of about 3.5 n . In Fig. [12] we display the result of our calculations 
for a collision of Au + Au at 2 AGeV. Comparing the left and the right hand part one 
recognizes the larger effect of the assumed pole mass shifts of the mesons in contrast to 
the C + C case shown in Figs. [8] and [91 Furthermore, it is clearly seen that the amount of 
di-electrons coming from p mesons is relatively larger than for light systems. This can be 
understood by the fact that during the longer collision time p mesons can rapidly decay 
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Figure 10: Comparison of the HADES di-electron spectrum with a calculation where the 
uj (p) cross section is scaled down by a factor 0.2 (0.8). 
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Figure 11: Transverse momentum spectra for three mass bins. Data source: [6. 
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Figure 12: Di-electron spectra for Au(2 AGeV) + Au. Left panel: In-medium spectral 
function is used and various contributions to the total spectrum are shown. Right panel: 
Vacuum spectral function (solid line) in comparison to the spectrum obtained with the 
in-medium spectral function (dashed violett line). 



and regenerate 

The effects discussed with respect to Fig. [H] are more clearly seen for a larger collision 
system. Figure [12] shows the dramatic change of the di-electron spectrum emitted from 
p and io mesons in central collisions Au + Au. The vacuum spectral function of the p 
meson (right panel) still shows the peak near the pole mass despite the m~ 3 dependence. 
In case of medium modification the shape differs strongly from the vacuum one. 

Repeating the same analysis as in Fig. [7] for central Au + Au collisions system we find 
that di-electrons from the dense region (dot-dashed lines) have low masses around 600 
MeV and contribute roughly 10% to the total u yield (see Fig. [T^l) . There is a remarkable 
difference between the outcome of the standard life time expression and Eq. (JT3J). 
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Figure 13: The same as in Fig. but for central Au + Au collisions. 
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Figure 14: As in Fig. [?] but for the heavy system Au+Au. 



4 Results at 1 AGeV 



Recently new HADES data are available also for C + C collisions at a bombarding energy 
of 1 GeV per nucleon [9]. The excess energy is about 450 MeV, thus only low energy 
tails of the p and u mesons play a role. Nevertheless the p mesons contribute essentially 
to the di-electron spectrum above an invariant mass of 500 MeV since other sources are 
even much smaller, see Figs. IT5landlT6l 

We obtain a reasonable agreement with the measured HADES [9] (Fig. IT51) and DLS 
[6J (Fig. [TBI data, but the data at an invariant mass around 400 MeV are underestimated. 
The shoulder in the data could only be explained by a higher contribution of di-electrons 
coming from the Dalitz decay of the rj mesons. We calculate a production cross section 
a v = 1.8 mb which has to be compared to the value a n o = 450 mb. The experimental 
values measured by [67], = 1.5 ± 0.4 mb and a n o = 287 ± 21 mb, give a ratio which 
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Figure 15: Di-electron invariant mass spectrum for C(l AGeV) + C calculated with 
in-medium spectral functions. Left panel: Spectrum in full phase space. Right panel: 
Comparison with HADES data J^j. 
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Figure 16: Comparison of the di-electron spectrum (line codes as in Fig. [73]) calculated 
with the matter spectral functions with the the measurements of the DLS collaboration 
The filter described on the DLS web page [7(tf was used. 



is a factor of 1.3 larger than our calculated value. Such an small increase of the rj yield 
could hardly improve the total di-electron spectra and explain the shoulder at 400 MeV. 

Figure fT7l exhibits the transverse momentum spectra for three invariant mass bins. A 
good agreement with available data can be stated. 

Remarkable is the following scaling property. Comparing the normalized di-electron 
spectra for the reaction n + p at 1.25 GeV and C + C at 1 GeV (see Fig. [TSJ) one 
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Figure 17: Transverse momentum spectra for three mass bins at 1 AGeV bombarding 
energy. Data source: [71^. 



recognizes that the spectra agree to a large extent despite of the different systems and 
beam energies. 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

M [GeV] M [GeV] 



Figure 18: Comparison of invariant mass spectra for n + p at 1.25 GeV (left panel, line 
codes as in right panel) and C + C at 1 AGeV (right panel). 



5 Summary 

In summary we have considered the propagation of broad resonances within a kinetic 
theory (transport) approach to heavy-ion collisions. Vector mesons are described by 
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spectral functions and these are evolved in space and time by a test-particle method. The 
motivation for this work is a new generation of data on di-electrons. The corresponding 
experiments are aimed at seeking for imprints of chiral symmetry restoration as particular 
aspect of in-medium modifications of hadrons. This lets us focus on the treatment of p 
and u mesons. The wildly wide-spread predictions call for an experimental clarification, 
but still heavy-ion data need often the comparison with models to extract the wanted 
information from data. 

We have utilized here the transport equations from Ref. which are approxima- 
tions of the much more involved Kadanoff-Baym equations [72]. Compared to an ap- 
proach wherein the spectral function is frozen in after creation the present framework 
let the spectral functions evolve towards the vacuum spectral functions. Therefore, the 
in-medium modifications are washed out, in particular, for the u meson. In contrast 
to earlier expectations the u peak does not suffer a significant modification, even when 
assuming a strong hypothetical shift of the peak position. Within the employed frame- 
work, medium modifications of p and to mesons are hardly seen in the di-electron spectra 
of small collision systems, even when using fairly strong and schematic assumptions for 
them. Only heavy collision systems seem to allow still to identify the wanted medium 
modifications. The elementary channels which contribute to the overall yields need better 
control to arrive at firm conclusions on interesting many-body effects. 
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